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ABSTRACT 

We reanalyse the measurements of parallax, proper motion , and line-qf-sight velocity 
for 18 masers in high mass star-forming regions presented bv lReid et alj (2009). We use 
a likelihood analysis to investigate the distance of the Sun from the Galactic centre, 
Ro, the rotational speed of the local standard of rest, vq, and the peculiar velocity of 
the Sun, Vq, for various models of the rotation curve, and models which allow for a 
typical peculiar motion of the high mass star-forming regions. 

We find that these data are best fit by models with non-standard values for v Q 
or a net peculiar motion of the high mass star-forming regions. We argue that a 
correction to v Q is much more likely, and that these data support the conclusion 
of IBinnevI (|2009l ) that Vq should be revised upwards from 5.2 km s -1 to llkms -1 . 
We find that the values of Ro and vo that we determine are heavily dependent on the 
model we use for the rotation curve, with model-dependent estimates of Rq ranging 
from 6.7 ± 0.5 kpc to 8.9 ± 0.9 kpc, and those of v$ ranging from 200 ± 20kms _ 
to 279 ±33kms" 1 . We argue that these data cannot be thought of as implying any 
particular values of Rq or vq. However, we find that vq/Rq is better constrained, lying 
in the range 29.9 — 31.6 km s~ kpc~ for all models but one. 
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1 INTRODUCTION 

The fundamental parameters that define the Solar position 
and velocity within the Galaxy remain uncertain to a re- 
markable degree. The major remaining uncertainty is that 
in the distance from the Sun to the Galactic centre, Ro- 
The most recen t results from studies of stellar orbits i n the 
Galactic centre (|Ghez et al.ll2008l; iGillessen et al.ll2009h give 
values of i? = (8.4 ± 0.4) kpc and R = (8.33 ± 0.35) kpc 
respectively. Th ese ca n be co mpared to the earlier estimate 
in the review bv lReidl l| 19931) of R = (8.0 ± 0.5) kpc. 

The total velocity of the Sun about the Galactic cen- 
tre is the sum of the velocity of the local standard of rest 
(LSR), vo, and the peculiar motion of the Sun with respect 
to the LSR in the same direction, V©. This total velocity 
can be determined using the apparent proper motion of Sgr 
A*, /Ua», since it is expected to be moving with a pecu- 
liar velocity less than ~ lkms -1 at the Galactic Centre. 
This constraint on the peculiar motion of Sgr A* is jus- 
tified by the observation that the velocity of Sgr A* per- 
pendicular to the pl ane is consistent with zer o, with uncer- 
tainties ~ lkms -1 (|Reid fc Brunthalerl \2004 ) . and because 
this motion is thought to be due to stochastic forces from 
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discrete interactions with individual stars, so the velocity 
in the plane should be similar to tha t perpendicular to it 
(IChatteriee. Herncmist. fc Loebl 12002t h iReid fc Brunthalerl 
(|2004h found lla* = (6.379 ± 0.024) mas yr" 1 , which cor- 
responds to (v + V Q )/Ro = (30.2 ± 0.2) km s " 1 kpc" 1 , or 
a velocity about the Galactic centre (using the iGhez et al.l 
result for Ro) of vo + Vq = (252 ± ll)kms _1 , where by 
far the dominant uncertainty comes from the value of Rp. 
Anal ysis of the GD-1 stellar stream l|Koposov. Rix. fc Hogd 
2009) has recently been used to suggest a significantly lower 
value, vo = 221^20 kms -1 . 

The velocity of the Sun with respect to the local stan- 
dard of rest, V0, is often assumed to be well known and con- 
strained to within ~ 0.5 kms" 1 because o f the analysis of the 
dynam ics of nearby stars conducted by iDehnen fc Binnevl 
henceforth DB98), and again more recen tly, us- 



ing identical techniques, by I Aumer fc Binnevl (|2009l ). DB98 
found 



v = (U Q ,V Q ,W G ) 

= (10.00 ± 0.36, 5.25 ± 0.62, 7.17 ± 0.38) kms" 1 , 



(1) 



which has been widely accepted and used. However IBinnevI 
|2009l . henceforth B09) suggests that the value for V© de- 
termined in these papers may be an underestimate by 
~ 6kms~ 1 . This is because the analysis by DB98 uses 
Stromberg's equation, which is derived under the assump- 
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tion that the Galactic potential is axisymmetric, and extrap- 
olates to zero velocity dispersion. In practice, the Galactic 
potential is not axisymmetric, and the smaller the velocity 
dispersion of a population, the more it is affected by depar- 
tures from axisymmetry. B09 uses a more global approach, 
ensuring that more emphasis is placed on stellar populations 
with high velocity dispersions, which one would expect to be 
less affected by the non-axisymmetry of the Galactic poten- 
tial. 

iReid et al.l (|2009T ) brought together observations of 
masers seen in high mass star-forming regions (hmsfrs) in 
the Milky Way. They used simple statistical tools in an effort 
to determine the values of Ro and vo, using the DB98 value 
of vq, initially under the assumption that the hmsfrs were 
moving on circular orbits in a flat rotation curve, for which 
they offer best-fitting parameters Ro — (8.24±0.55) kpc and 
vo — (265 ± 26) kms -1 , but with a high \ 2 value. 

In addition they considered a model in which the 
hmsfrs were moving with a characteristic velocity with re- 
spect to their local circular velocity. They found that they 
achieved a significantly improved fit to their data using this 
model with a peculiar velocity ~ 15 kms -1 in the opposite 
direction to rotation. This fit yielded Ro — (8.40 ±0.36) kpc 
and vo = (254± 16) kms -1 . They briefly considered a model 
in which the value of v© was allowed to vary, finding an ac- 
ceptable fit with v© = (9, 20, 10) kms -1 . However, believing 
the DB98 result to be "well determined" they did not pursue 
the matter further. 

In thi s pap er we re-examine the data described in 
IReid et Ifll l|2009T ). and conduct a likelihood analysis for var- 
ious models of the velocity distribution of the maser sources. 
This enables us to exploit these data more thoroughly. We 
also consider the implications of the results from B09 on the 
interpretation of these data. In Section [2] we explain what 
the data consist of, describe our models and our statistical 
technique; in Section we give the raw results and exam- 
ine their significance. We discuss the implications of these 
results in Section U 



2 METHODS 
2.1 The data 

These data, as given in Table 1 of lReid et all (|2009h . consist 
of measurements for 18 masers of: Galactic coordinates (h, 
bi), which we can assume to be exact; parallaxes 7r;; proper 
motions and jJ, y ,i\ line-of-sight velocities vlsr.i- For each 
quantity Reid et al. give an error and we assume that this 
error together with the measured value of the quantity de- 
fines a Gaussian probability distribution for the true value of 
the quantity. For the proper motions [t, x and [i y we shall use 
the more conventional notation \i a and The line-of-sight 
velocity «lsr is relative to an obsolete estimate of the LSR. 
Fortunately the underlying heliocentric line-of-sight veloc- 
ity v r can be recovere d from fLSR with out impact on the 
associated uncertainty (|Reid et alj|2009l . Appendix). 



direction of rotation (v ) 




X (kpc) 

Figure 1. Diagram showing a pole-on view of the Galaxy il- 
lustrating the various peculiar velocities we consider (e.g. equa- 
tion This therefore only shows the in plane components. The 
Sun (represented by a solid circle) is placed at (0, 8) kpc, and 
components of its velocity vq are indicated. Three points (empty 
squares) are plotted to represent masers, and the directions of 
the components of any systematic offset from circular velocity in 
each case (vsfr) are shown, as is the peculiar velocity due to any 
bias in the observed radial velocity, %. This diagram is purely 
illustrative and does not show any real data. 

2.2 Our models 

The motion of both the Sun and the masers is dominated 
by circular motion around the Galactic centre with velocity 

-v c (R)e 4 , (2) 

where v c (R) > and the minus sign reflects the fact that 
the Galaxy rotates clockwise in our coordinate system. We 
explore three forms for v c (R) 

• A flat rotation curve, v c (R) = Vq, with vo being a (pos- 
itive) free parameter. 

• A power-law rotation curve v c (R) — Vo{R/ Ro) a , with 
Wo, -Ro and a being free parameters. 

• A rotation curve corresponding t o that given by the 
Galac tic potential Model I in §2.7 of iBinnev fc Tremaind 
(2008, henceforth GDII), linearly scaled to variable values 
of Ro and va- 
in each case, we take the parameters of v c (R) to have uni- 
form prior probability distributions. This choice ensures that 
we determine what these data tell us, without prejudice. 

The velocity of a maser can be expected to differ from 
the local circular speed. We separate this difference into a 
random comp onent and ( i n som e cases) two systematic com- 
ponents. Like lReid et al.l l|2009t ). we consider the possibility 
that the velocity of hmsfrs has a systematic offset from the 
circular velocity, vsfr- We also consider the possibility that 
the expansion of the shell within which a maser occurs will 
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displace the maser's velocity from that of the star, and if we 
are biased towards seeing masers on the near side of that 
shell, we will see a bias in radial velocity. We represent this 
as an average peculiar motion w m e«, where e* is the unit 
vector from the Sun to the maser. We assume that the ran- 
dom component of the maser velocities has a Gaussian dis- 
tribution of dispersion A„. In total we take the probability 
distribution of the velocities of masers to be 



p(v) d v 



where 



d 3 v 



(2ttA2) 3 / 2 



exp 



2A2 



v = -v c (R)e^ - vsfr + v m e*. 



(3) 



(4) 



Notice that positive values of «</,,sfr. imply that the masing 
stars lead Galactic rotation, and that motion of the masers 
towards us will be reflected in negative values of v m . In most 
cases we fix both vsfr = and v m = 0. The various veloci- 
ties are illustrated in Fig. [T] 

We sometimes take the Sun's motion with respect to 
the LSR, V0, to be specified, and sometimes we fit it to the 
data. 

For given values of h, bi, the heliocentric distance Si, 
and Ro, we deduce the probability distributions in proper 
motion and line-of-sight velocity given the above probability 
distribution of the velocity (eq.0. Since our focus is on the 
Galaxy's parameters, we marginalise over the Si. We do this 
under the assumption that the probability distribution of 
the Si is that implied by the assumption of Gaussian errors 
in the parallaxes. 



2.3 Statistical analysis 

2. 3. 1 Likelihood function 

In order to determine the best-fitting model parameters, we 
maximise the likelihood function. This is 



C(0) oc 



p(data | 0), 



(•») 



where p(data | 0) is the conditional probability of the ith 
observation, given the model with parameters 



= (v ,Ro, a, v Q , A„, Si,...). 



(6) 



For each distance Si, the model provides a multivariate 
Gaussian probability distribution in v r ,i, fJ, a ,% and fj,gj. 
The data also provide a multivariate Gaussian distribu- 
tion for the underlying values of these variables. We obtain 
p(data| 0) by integrating the product of these two Gaussian 
distributions over all three variables. 

The likelihood function defines the a posteriori prob- 
ability distribution o f the model parameter s 0. We use a 
Metropolis algorithm ([Metropolis et aU1953T ) to identify the 
peak in this probability distribution, and to characterise its 
width around the most probable model. The Metroplis algo- 
rithm is a Markov Chain Monte Carlo method for drawing a 
representative sample from a probability distribution, such 
as the likelihood function. We start with some choice for the 
parameters 0, and calculate the associated likelihood C{0). 
We then 

(i) choose a trial parameter set 0' by moving from in 



all directions in parameter space, by an amount chosen at 
random (from a Gaussian distribution); 

(ii) determine C{0'); 

(iii) choose a random variable r from a uniform distribu- 
tion in the range [0,1]; 

(iv) if C(0') j C(0) > r, accept the trial parameter set, and 
set = 0'. Otherwise do not accept it. 

(v) Return to step (i). 

The first few values of are ignored as "burn-in", which 
helps to remove the dependence on the initial value of 0. We 
repeat the procedure until the chain of values constitutes 
a fair sample of the probability distribution - we establish 
that the burn-in period is sufficiently lo ng by comparing 
chain s that have different starting (e.g. iGelman fc RubirJ 

Eli). 



2.3.2 Bayesian evidence 

In addition to comparing models that differ only in the val- 
ues taken by a given set of parameters, we have to assess the 
value of adding a parameter to a model. Adding a parameter 
is guaranteed to increase the maximum likelihood achievable 
for given data, but is that increase statistically significant 
or just the consequence of an enhanced ability to fit noise? 
The Bayesian methodology for making such ass e ssmen ts is 
now well established and described by iHeaven 1 11200911 . for 
example. 

One calculates the "evidence" for the model under dif- 
ferent priors for whatever parameters are either fixed a priori 
or varied. The evidence is the total probability of the data 
after integrating over all parameters: 

p(data| Model) = j d n 0p(data|0, Model) p(0jModel). (7) 

Here p(0|Model) is the prior on the parameters. For a pa- 
rameter 9i that is varied, we take the prior to be uniform 
within a range of width A6i that is large enough to encom- 
pass any plausible value of 6i (so p(0i\ Model) A9i = d9i/ A9i 
over this range). The priors on parameters that are always 
varied play little role because we are interested the ratio of 
the evidence when a parameter 9 n is fixed to when it is var- 
ied. When On is fixed, its prior is a delta function at the 
chosen value, so the ratio of the evidences when 9 n is fixed 
to when it is varied is 



p(data| Model, 0„ fixed) _ / d"- 1 6>p(data|0, Model) 



p(data| Model) 



/ d™0p(data|0, Model) 



A0„,(8) 



where upper integral excludes 8 n and the lower one includes 
it. 

By Bayes theorem the ratio of the a posteriori proba- 
bilities of the model when 8 n is fixed to when it is varied is 
the ratio of the corresponding evidences times the ratios of 
the prior on 9 n being fixed to that on it being variable. We 
take this second factor to be unity, so the a posteriori prob- 
abilities of the models is simply the ratio of the evidences. 

The integration over the space of parameters that is re- 
quired to calculate the evidence is exceedingly costly if done 
by brute force. We approximate it by assuming that in the 
vicinity of its peak, p(data|0, model) can be approximated 
by a Gaussian p oc exp(— T ■ K ■ 0), where if is a matrix 
whose eigenvectors and eigenvalues can be estimated from 
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Figure 2. Plot showing contours of the Likelihood function 
(marginalised over A„) for a model with the GDII rotation curve, 
and the DB98 vg . There is clearly a strong correlation between 
the values found for vq and Ho . Contours are drawn at likelihood 
differing from the maximum likelihood AC = 0.25, 0.5, 1, 2. 



the output of the Metropolis algorithm. With this approxi- 
mation, the integral over parameters becomes analytic. 

The Metropolis algorithm yields a set of points in pa- 
rameter space which sample p(data| Model). Principal com- 
ponent analysis of this sample yields the eigenvectors and 
eigenvalues of K. 



3 RESULTS 

We first investigate models in which the masers are at rest 
with respect to their host stars (i.e., Vm = 0). 

Table Q] gives the peak likelihoods of our models, as 
well as the best-fitting parameters and the corresponding 
uncertainties. The best-fitting value of vo varies in the range 
(200 - 279) kms -1 , and that of R between 6.7 and 8.9kpc. 
As Fig. [2] illustrates, Vo and Ro are strongly correlated with 
the result that vo/Ro is confined to the relatively narrow 
range 29.9 — 31.6 kms -1 kpc -1 . 

When vq is a free parameter, the best-fitting values of 
Uq and Wq are close to the DB98 values, and essentially 
independent of the form of the rotation curve, whereas Vq 
varies in the range 16.5 — 19.5 kms" 1 . Similarly, when vq is 
fixed at the DB98 value and vsfr is a free parameter, only 
f0,SFR takes a value that is far removed from that which 
one would naively expect - it moves in the range —11.0 to 
— 14.8 kms -1 . Thus the data suggest either that the Sun 
is circulating significantly faster than the circular speed, or 
that the hmsfr have an appreciable rotational lag. 

The fits are almost perfectly degenerate between vq 
and vsfr: they constrain only the difference between these 
velocities. However, we shall argue in Section [4] that signif- 
icantly non-zero values of vsfr are physically implausible, 



so here we focus on what can be inferred about vq given 
vsfr. = 0. In view of the degeneracy between vq and vsfr 
we do not report results obtained when both vq and vsfr 
are varied. 

When vsfr = 0, the peak likelihood is higher when vq 
is fixed to the B09 value than when it is fixed to the DB98 
value. To assess the significance of this increase in likelihood, 
we calculate the ratio of the evidences for the two models as 
described in Section \2. 3. 21 

f 2 x 10 -4 for a = 

1 x 10 -4 a variable (9) 
[ 4 x 10 -4 GDII. 



p(DB98|data) 
p(B09|data) 



Thus regardless of the adopted rotation curve, the data 
strongly favour upward revision of Vq from 5.2 km s -1 to 
11 kms -1 . 

When vq is a free parameter, the peak likelihood of the 
B09 model is surpassed at yet larger values of Vq. To de- 
termine whether the increase in likelihood that occurs when 
vq is set free from the B09 value, we again calculate the rel- 
evant ratio of the evidences. Since the two models now differ 
in whether v is fixed or free, the priors on the components 
of vq now become relevant (cf eq. [8]). We have adopted 
AU Q = AVq = AW Q = 100 kms -1 . With these values we 
have 



p(B09jdata) 
p(vq free|data) 



5 for a — 
0.6 a variable 
40 GDII. 



(10) 



Therefore under these assumptions the increase in likelihood 
attained on setting vq free from the B09 value is not statis- 
tically significant. The ratio of evidences for the case when 
vq is set to the DB98 value or is set free is 




p(vq freejdata 



for a = 
a variable 
GDII. 



(11) 



Therefore, even with a generous choice of AUq, etc., the 
data reject the possibility that the Sun has the DB98 value 
of vq. Reducing the widths AUq, etc., of the priors on vq 
would strengthen the case against the DB98 value of Vq. 

The choice of A Uq ~ 100 kms -1 is reasonably gen- 
erous, and it is sensible to ask what value of AUq, etc., 
would bring p(B09|data)/p(vQ free|data) down to unity. In 
the case where the GDII rotation curve is used, it would 
require a reduction to a value AUq ~ 30 kms -1 - this is 
smaller than is reasonable given that the value for Vq found 
when it is allowed to vary is already ~ 15 kms -1 greater 
than the canonical DB98 value. 

It is worth noting that the case for setting vq free from 
either the DB98 value of the B09 value is weakest when the 
best-motivated rotation curve is adopted - the GDII curve. 

The lowest values of both vo and Ro are found for the 
models with the DB98 value of v , and vsfr set to zero. 
Both vo and Ro take larger values for models with the B09 
value of vq, and larger still for models with either vq or 
vsfr allowed to vary. The lowest values of vo and Ro for the 
models in which the peculiar velocities are allowed to vary 
are 232 kms - and 7.7 kpc respectively. 

Fig. 13 suggests a reason why the value of vq or vsfr 
in the model has such a large impact on the best fitting 
values of vo and Ro- It shows the residual velocities of the 
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Table 1. Log likelihoods for best-fitting models with fiat rotation curves (a = 0), power law rotation curves (a 7^ 0) and rotation curves 
taken from GDII. The first 9 likelihoods are for models in which the sources are assumed to have no typical velocity - for the first 3 the 
value of vq is that found in DB98; the next 3 are for models in which vq is assumed to be that suggested by B09; and the next 3 are for 
models in which vq is allowed to vary freely; the final 3 are for cases in which vq is taken to be the DB98 value, and the mean source 
peculiar motion with respect to their local standard of rest, vgpR, is allowed to vary freely. 
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Ro 
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a 


(kms" 1 ) 


(kms" 1 ) 


Wq 
(kms" 1 ) 


(kms^ 1 ) 


v /Ro 
(km s — x kpc _ 


') 


log(£) 


DB98 


onn _l_ on 
2U0 ± 2U 

209 ± 26 

218 ± 20 


6.7 ± 0.5 
6.9 ±0.7 
7.1 ±0.5 




0.10 ±0.16 

f^FlTT 

KjlJlL 


10.0 
10.0 
inn 


5.2 
5.2 

•J.Z 


7.2 
7.2 
t .z 


10.0 ± 1.3 
10.0 ± 1.3 
9.3 ±1.2 


30.1 ± 1.7 
30.4 ± 1.7 
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26.8 
27.1 

Qfl 1 
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B09 
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7.0 ±0.5 
7.4 ±0.6 
7.6 ±0.5 




0.15 ±0.13 
GDII 


10.0 
10.0 
10.0 


11.0 

11.0 
11.0 


7.2 
7.2 
7.2 


8.1 ± 1.1 
8.0 ±1.1 
7.6 ±1.0 


30.5 ± 1.5 
30.8 ± 1.5 
31.1 ± 1.5 




35.7 
36.5 
38.9 


m 

C> 

> £ 


232 ± 24 
258 ± 32 
246 ± 24 


7.7 ±0.6 
8.6 ±0.9 
8.1 ±0.6 




0.27 ±0.13 
GDII 


8.1 ±2.9 
8.1 ± 2.8 
8.3 ±2.8 


18.6 ±2.4 
19.5 ± 2.5 
16.5 ±2.4 


9.7 ± 2.0 
10.1 ±2.0 
9.9 ±2.0 


7.1 ± 1.0 
7.1 ± 1.0 
7.0 ±1.0 


30.0 ± 1.8 
29.9 ± 1.7 
30.3 ± 1.8 




42.5 
45.4 
43.5 




-co 


Ro 


a 


V R,SFR 


"0,SFR 


«z,SFR 


A,. 


vo/Rq 




log(£) 


a 

b CD 

CO CD 

> a 


241 ± 24 
279 ± 33 
259 ± 23 


7.7 ±0.6 
8.9 ±0.9 
8.2 ±0.6 




0.25 ±0.12 
GDII 


3.3 ±2.8 
2.0 ±2.8 
2.5 ± 2.8 


-12.9 ±2.4 
-14.8 ±2.5 
-11.0 ±2.4 


2.5 ±2.0 
3.0 ±2.0 
2.8 ±2.0 


7.0 ±1.0 
6.6 ±1.0 

7.1 ± 1.0 


31.1 ± 1.7 

31.5 ± 1.6 

31.6 ± 1.7 




42.9 
45.5 
43.8 


t>o( km 


s" 1 ) i? (kpc) 


Uq 


^0 


Wq 


% 


A„ Do/i?o(kms kpc" 




log(£) 


233 ±20 7.3 ±0.5 10.0 5.2 7.2 -8.1 ±2.7 8.0 ±1.1 
247 ±19 7.8 ±0.4 10.0 11.0 7.2 -6.2 ± 2.4 6.8 ± 0.9 
259 ± 25 8.2 ±0.5 10.7 ±3.3 15.0 ± 2.7 10.0 ± 2.0 -4.8 ± 2.9 6.5 ± 1.0 


31.6 ± 1.5 
31.4 ± 1.4 
31.4 ± 1.9 




36.2 
43.1 
46.0 


tio(km 


S 1 ) iJo(kpc) ffl,SFR v 


A.SFR 


»z,SFR 




A„ iio/ft(kms kpc" 




log(£) 


271 ±25 8.2 ±0.6 2.0 ±3.1 -9.7 ± 2.8 - 


2.5 ±2.2 -5. 


8 ±3.3 6.4 ±1.1 


32.9 ±2.0 




46.0 



Table 2. Similar to TableQ"] This shows the best-fitting parameters (and corresponding log likelihoods) for models with a GDII rotation 
curve, and with the value i> m allowed to vary. 

objects after the expected velocity is subtracted (ignoring 
the uncertainty in parallax). The objects that provide the 
strongest indication that some change in peculiar motion is 
required (indicated by open squares) are all relatively close 
to the Sun. Consequently, changes in i?o and vq have a rel- 
atively small "lever arm" with which to bring the model 
closer to the data, with rather large changes required in or- 
der to have any significant impact on the expected values of 
the observables for those objects. The introduction of a new 
peculiar velocity (either of the sources or the Sun) allows 
the model to compensate for the significant discrepancy be- 
tween the velocity of these objects and those expected from 
the rotation curve, without producing a major effect on the 
best-fitting values of Rq and vo- 



3.1 Maser line-of-sight velocities 

Finally we consider models in which we allow for a system- 
atic difference between the measured line-of-sight velocity of 
a maser and the line-of-sight velocity of the exciting star - 
such a difference would arise if masing occurred on the near 
side of an expanding shell around the star. We have done 
this for all the rotation curves described in this paper, but 
the results are sufficiently similar to one another that we 
only present (in Table [2} the results for the GDII rotation 
curve. 



We again compare the different models by calculating 
the ratios of their evidences, using the prior for v m uniform 
in a range of width Av m = 20kms -1 . We find 



p{vm = 1 data) 
p(vm free|data) 



0.01 for DB98 

0.06 B09 . . 

0.3 v Q free ( ' 

0.3 vsfr free 



so in each case the data support adding the extra free pa- 
rameter. 

If we accept this extra parameter in our models, we have 
to reconsider the evidence for a change in vq (or vsfr). The 
relevant ratios of evidences are: 

p(DB98,u m free|data) _ Q QQ1 

0.3, (14) 



p(B09, v n 


free|data) 


p(DB98,-u 


t! free|data) 


P(V , Vm 


free|data) 


p(B09,u m 


free|data) 



, i , , - 220. (15) 
p(v0,w m free|data) 

So the B09 value of v Q is still favoured. Moreover, when 
vq is taken to be free, the likelihood still peaks at an even 
larger value of V than that of B09. 

We note that the sources providing the strongest evi- 
dence that a non-zero value of v m is needed are those ringed 
in the the bottom right plot of Fig. [3] all associated with the 
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Figure 3. Residual velocities left after the best-fitting model velocities are subtracted, for models with GDII rotation curves and vq 
taking the DB98 value (top-left); the B09 value (top-right); or taken to be a free parameter (bottom- left). We also plot the residual 
velocities left after subtracting the best fitting model velocities for a model with the B09 vq value, and « m taken to be a free parameter 
(bottom-right). Open squares correspond to the objects for which the likelihood improves most significantly between the best-fitting 
DB98 model and the best-fitting free Solar velocity one (A(log£i) > 1), and crosses to all other sources. The Sun is represented by a 
solid circle. The solid lined ellipses around each point are 1 a measurement error ellipses, where the velocity error perpendicular to the 
line of sight in each case is due to a combination of the uncertainty in the proper motion and in the parallax. The uncertainty in the 
parallax also means that the residual velocities shown here are not an ideal illustration of the difference between the model and the data, 
because the peculiar motions can only be given for a chosen position (in this case the position corresponding to the quoted parallax). 
However this does provide a useful guide, especially for sources relatively close to the Sun, which have small position uncertainties. The 
dashed ellipse in the bottom-right plot is drawn around the three objects for which the likelihood improves most significantly when ri m 
is taken to be a free parameter (A(log£i) > f for B09 Vq) 



Perseus spiral arm (Re id et al.ll2009h . If we exclude these 
sources, the best fitting value of v m is ~ — 3kms _1 for any 
assumption about vq we consider. This is approximately 
within the uncertainty on v m , and if we considered only 
this subset of the data in equation ([12} , the evidence would 



not support adding the extra free parameter v m (though it 
should be noted that it is still a radial velocity towards the 
Sun). It is, therefore, possible that what we have modelled 
as an offset in the radial velocities of all observations is actu- 
ally primarily due to a large peculiar velocity of the objects 
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in the Perseus arm, directed approximately in the direction 
of the Sun. 



4 DISCUSSION 

These results, and in particular the large variation in vo 
and Rq depending on the other model parameters, makes it 
impossible to constrain tightly either vo or Ro from these 
data - the smallest best-fitting values of these parameters 
are, respectively, 40 per cent and 33 per cent smaller than 
the largest best-fitting values. The choice of rotation curve 
has a significant impact, and we do not investigate all pos- 
sible rotation curves, ft is worth noting that our two most 
favoured models in Table [1] (according to Bayesian evidence, 
so somewhat dependent on our choice of priors) have best- 
fitting values of vo that differ by 30 kms -1 and of Ro that 
differ by f kpc from one another (the models are those with 
B09 v© + GD rotation curve and free vq + power-law ro- 
tation curve, respectively). 

However, for all but one of our models the best fitting 
vo/Ro lie in the narrow range (29.8 — 31.5) kms -1 kpc -1 , 
with typical uncertainties ~ 1.5 kms -1 kpc -1 . This corre- 
sponds to best fitting (no + Vq)/Ro in the range (30.9 — 
32.5) kms -1 kpc -1 (again, depending on the model), with 
similar uncertainties. These values are slightly larger than 
the value of (v + V Q )/R = 30.2 ± 0.2 kms -1 kpc -1 foun d 
from the proper motion of Sgr A* (|Reid fc BrunthalerT 2004). 
but consistent to within (less than) twice the uncertainties 
on our values. 

If we allow for a bias in the maser radial velocities with 
respect to the Sun by allowing the parameter v m to take 
a non-zero value in equation Q, this does improve the fit. 
However we have seen that this is primarily due to a group 
of maser sources in the Perseus spiral arm, so the cause may 
be a large peculiar motion in that part of the Galaxy. Even 
if the bias in measured radial velocity is real, it does not 
affect the result: (i) that in different models these data are 
fit best by very different values of vo and Ro; (ii) that the 
DB98 vq is rejected by these data, and (iii) that the B09 
V0 is favoured o ver m aking v© a free parameter. 

iReid et all (|2009h focused on the idea that these data 
show that the hmsfrs are orbiting the Galaxy with a ve- 
locity that lags the circular speed by — u^.sfr ~ 15 kms -1 . 
Our analysis of similar models shows that the data are fit 
better by a slightly smaller offset of (11 — 14.8) kms -1 to 
the circular speed. We have also shown that an equiva- 
lent fit to the data can be obtained by instead increasing 
to ~ 17 kms - the amount Vq by which the Sun is as- 
sumed to circulate faster than the circular speed. In fact, 
the sources that provide the strongest statistical support for 
iV>,sfr ~ —15 kms -1 are found near the Sun (Fig. [3]), so a 
change in vq has a very similar effect to a change in vsfr. 
More maser data at different Galactic azimuths could break 
this degeneracy, as well as reducing the correlation between 
the values of vo and Ro seen in Fig . [2] 

If the proposal of IReid et aU ([2009|) that i^,sfr — 
— 15 kms -1 were correct, the hmsfrs would all have to be 
close to apocentre. Since the lifetime of an hmsfr is short 
compared to a typical epicycle period and spiral structure 
must play a significant role in their formation, it is not in- 
herently implausible that the maser stars are all close to 



apocentre. However, two arguments make it unlikely that 
the maser stars are on orbits as eccentric as is implied by a 
15 kms -1 offset from circular motion at apocentre. 

First, non-axisymmetric structure in the Galaxy's po- 
tential modulate s the tangential velocity of gas by only 
~ 7kms -1 (e.g. iBinnev fc Merrifield1ll998l . §9.2.3), so the 
hmsfrs would be moving with a peculiar velocity signifi- 
cantly higher than the gas from which they formed. Sec- 
ond, by GDII (equation 3.100), a quarter of an epicy- 
cle period later the U velocities of these stars would be 
(2f2//t)15 kms -1 ~ 22kms -1 , where O and re are the cir- 
cular and radial frequencies at the star's location. Hence 
the radial velocity dispersion of a population of such stars 
would be at least ~ 22/^/2 ~ 15 kms -1 . The velocity dis- 
persion of stars observed locally increases with age on ac- 
count of heating by irregularities in the Galaxy's gravita- 
tional field (e.g. GDII §10.4.1) and the velocity dispersion 
of the bluer stars in the Hipparcos catalogue is ~ 10 kms 
jAumer fc Binnevl 120091') rathe r than 15 kms -1 . Thus the 
conjecture of Reid et all (|2009T l not only requires the mas- 
ing stars to be confined to apocentre but also requires them 
to have more eccentric orbits than the generality of young 
stars. Recognising this problem, Reid et al. suggested that 
the orbits of these stars became more circular early in their 
lives. However, any random scattering process spreads stars 
more widely in phase space and therefore increase the mean 
eccentricity of the masing stars. Only a dissipative process 
could increase the phase-space density of these stars by mov- 
ing them to more circular orbits, and no such process is 
known. 

If the B09 value of Vq is correct, the lag to circular 
motion required to optimise the fit to the data is so small 
(~ 5 kms -1 ) that the above objections to the proposal that 
hmsfr systematically lag rotation become moot. However, 
the formalism of Bayesian inference says that when the B09 
value of Vq is accepted, there is no convincing evidence for 
a systematic lag of the hmsfr. 

Throughout this paper we have proceeded under the as- 
sumption that the velocity of the LSR is the same thing as 
the circular velocity at Ro- It is worth noting that this is 
not necessarily the case. The LSR is defined to be the veloc- 
ity of a closed orbit as it passes the current Solar position, 
which will only be a circular orbit if the Galactic potential is 
axisymmetric. Therefore the velocity of the LSR may be off- 
set from the circular velocity curve. However (much like the 
possibility of a systematic lag of the hmsfr) we should note 
that if we accept the B09 value of Vq, there is no convincing 
evidence for this offset. 

We conclude that these data provide a compelling case 
for revising Vq upward. Our Bayesian analysis (equation llO|l 
supports the revision upwards of Vq to 11 kms -1 suggested 
by B09 on the basis of modelling predominantly old disc 
stars with relatively large random velocities. It should, how- 
ever, be recognised that an even better fit to these data 
comes from somewhat larger upward revisions to Vq ~ 
16 - 20 kms -1 . 



5 CONCLUSIONS 

We have reanaly sed observations of hmsfrs reported in 
IReid et al.l (|2009f ) using a maximum-likelihood approach to 
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exploit fully the data, in an effort to determine the val- 
ues of Ro, the Galactocentric radius of the Sun, and vo, 
the circular velocity of the LSR. We have found that the 
best-fitting values, considered separately, are strongly de- 
pendent on the Galaxy model used to interpret the data, 
but that the ratio vq/Ro is consistently found to lie in the 
range 29.8 — 31.5 km s _1 kpc -1 . 

We have also used these data to explore the value of the 
Sun's peculiar velocity vq, in light of the recent argument 
of B09 that the canonical DB98 value is incorrect. We find 
that these data support the conclusion that Vq is signifi- 
cantly higher than the DB98 value. By a small but signifi- 
cant margin, the data prefer models with the B09 revision 
of Vq to ~ llkms -1 over models in which vq is left as a 
free parameter f equation 1 1 Op . The best- fitting models have 
Uq and Wq near to the DB98 values, and the even larger 
valu e of Vq ~ 16 - 20 km s" 1 . 

iReid et all (2009) suggested that hmsfrs significantly 
lag circular rotation. We have investigated the possibility 
that the hmsfrs have a typical peculiar velocity vspr and 
find that the models only constrain the velocity difference 
vsfr — vq. We have argued that models in which the hmsfrs 
have large peculiar velocities in the opposite direction to 
Galactic rotation are neither needed nor plausible. 

This work, in conjunction with B09, casts severe doubt 
on the accuracy of the widely used DB98 value for Vq. This 
must be a concern for anyone interested in the dynamics 
within the Milky Way because velocities are inevitably mea- 
sured with respect to the Sun. Both observations of more 
masers and the development of a detailed dynamical model 
of the Galaxy's spiral structure would contribute to estab- 
lishing more securely what the true value of Vq is. 



Ghez A. M., Salim S., Weinberg N. N., Lu J. R., Do T., 
Dunn J. K., Matthews K., Morris M. R., Yelda S., Becklin 
E. E., Kremenek T., Milosavljevic M., Naiman J., 2008, 
ApJ, 689, 1044 

Gillessen S., Eisenhauer F., Trippe S., Alexander T., Genzel 
R., Martins F., Ott T., 2009, ApJ, 692, 1075 

Heavens A., 2009, arXiv: 0906.0664 

Koposov S. E., Rix H.-W-, Hogg D. W., 2009, ApJ, sub- 
mitted, (arXiv: 0907.1085) 

Metropolis N., Rosenbluth A. W., Rosenbluth M. N., Teller 
A. H., Teller E., 1953, Journal of Chemical Physics, 21, 
10871092 

Reid M. J., 1993, AnnRA&A, 31, 345 
Reid M. J., Brunthaler A., 2004, ApJ, 616, 872 
Reid M. J., Menten K. M., Zheng X. W., Brunthaler A., 
Moscadelli L., Xu Y., Zhang B., Sato M., Honma M., Hi- 
rota T., Hachisuka K., Choi Y. K., Moellenbrock G. A., 
Bartkiewicz A., 2009, ApJ, 700, 137 
Rygl K. L. J., Brunthaler A., Reid M. J., Menten K. M., 
van Langevelde H. J., Xu Y., 2009, A&A, accepted 
(arXiv:0910.0150) 



ACKNOWLEDGMENTS 

We thank John Magorrian and the other members of the 
Oxford dynamics group for valuable discussions. PJM is sup- 
ported by a grant from the Science and Technology Facilities 
Council. 



NOTE ADDED AFTER ACCEPTANCE 

After this paper was accepted for publication, iRvel et al.l 
(2009) presented similar observations of maser sources. In- 
corporating these observations in our analysis does not ma- 
terially affect our conclusions. 
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